Stopping effects in U+U collisions with a beam energy of 520 MeV/nucleon 
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A Relativistic Transport Model (ART1.0) is applied to simulate the stopping effects in tip-tip 
and body-body U+U collisions, at a beam kinetic energy of 520 MeV/nucleon. Our simulation 
results have demonstrated that both central collisions of the two extreme orientations can achieve 
full stopping, and also form a bulk of hot, dense nuclear matter with a sufficiently large volume 
and long duration, due to the largely deformed uranium nuclei. The nucleon sideward flow in 
the tip-tip collisions is nearly 3 times larger than that in body-body ones at normalized impact 
parameter b/b max < 0.5, and that the body-body central collisions have a largest negative nucleon 
elliptic flow V2 = —12% in contrast to zero in tip-tip ones. Thus the extreme circumstance and the 
novel experimental observables in tip-tip and body-body collisions can provide a good condition and 
sensitive probe to study the nuclear EoS, respectively. The Cooling Storage Ring (CSR) External 
Target Facility (ETF) to be built at Lanzhou, China, delivering the uranium beam up to 520 
MeV/nucleon is expected to make significant contribution to explore the nuclear equation of state 
(EoS). 

PACS numbers: 24.10.Lx,25.75.Ld,25.75.Nq,24.85.+p 



I. INTRODUCTION 

In recent years, the ultra-relativistic high energy heavy 
ion collisions performed at SPS/CERN and RHIC/BNL 
(V s nn ~ 10 — 200 GeV) focus on high temperature 
and low baryon density region in nuclear matter phase 
diagram [l[ to search a new form of matter with par- 
tonic degree of freedom-the quark-gluon plasma (QGP) 
@j Hi EL H[ ■ However, no dramatic changes of experimen- 
tal observables, such as jet-quenching, elliptic flow and 
strangeness enhancement, have been observed yet [6[. On 
the other hand, the heavy ion collisions performed at the 
BEVALAC/LBNL and SIS/GSI 0,1 in last two decades 
were used to produce hot and compressed nuclear mat- 
ter to learn more about the nuclear equation of state 
(EoS) [HI, 03] at high baryon density and low temper- 
ature region of the phase diagram. Although we have 
made great efforts to study the nuclear EoS, theoreti- 
cally and experimentally, a solid conclusion can hardly 
be made. Then, it is still worthwhile to systematically 
study on the collision dynamics as well as the EoS observ- 
ables. Recently, for more understanding of the nuclear 
matter phase diagram and EoS at high net-baryon den- 
sity region, it is proposed to collide uranium on uranium 
target at External Target Facility (ETF) of Cooling Stor- 
age Ring (CSR) at Lanzhou, China with a beam kinetic 
energy of 520 MeV/nucleon. [nj. 

Uranium is the largest deformed stable nucleus, and 
has approximately an ellipsoid shape with the long and 
short semi- axis given by Ri = Rq(1 + 26/ 3) and R s — 
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FIG. 1: (Color online) (a) body-body collisions (b) tip-tip 
collisions 



Rq(1 — <5/3), respectively, where Rq = 7 fm is the effec- 
tive spherical radius and 5 = 0.27 is the deformation pa- 
rameter ||. Consequently, one has Ri/R s — 1.3. In our 
simulation, we consider two extreme orientations: the so- 
called tip-tip and body-body patterns with the long and 
short axes of two nuclei are aligned to the beam direc- 
tion, respectively [13], see Fig. [T] for illustration. The 
two types of orientations can be identified in random ori- 
entations of U+U collisions by making proper cutoffs in 
experimental data, such as the particle multiplicities, el- 
liptic flow and so on [l(J, [Tl[ • With the two extreme 
collision orientations, some novel stopping effects which 
are believed responsible for some significant experimental 
observables, such as particle production, collective mo- 
tion as well as attainable central densities, can be ob- 
tained. Due to the large deformation of the uranium 
nuclei [IJ, [HI , it is expected that the tip-tip collisions 
can form a higher densities nuclear matter with longer 
duration than in body-body or the spherical nuclei colli- 
sions, which is considered to be a powerful tool to study 
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the nuclear matter phase transition at high baryon den- 
sity , and the body-body central collisions may reveal 
a largest out-of-plane elliptic flow (negative V2) at high 
densities, which can be a sensitive probe to extract the 
early EoS of the hot, dense nuclear matter [f2|, H3|- The 
novel experimental observables can be effectively utilized 
to study the possible nuclear matter phase transition and 
the nuclear EoS [H G3, [H, QJ, [H, till EI GJ, HI . For 
comparing with tip-tip and body-body collisions, a type 
of gedanken "sphere-sphere" collisions without deforma- 
tions of uranium nuclei are also included in the simula- 
tion. 

The ART1.0 model [U H3] derived from Boltzmann- 
Uehling-Uhlenbeck (BUU) model [23[ has a better treat- 
ment of mean field and Pauli-Blocking effects [23| than 
cascade models (24{. The fragments production mech- 
anism and partonic degree of freedom are not present 
in the ART1.0 model. A soft EoS with compressibility 
coefficient K = 200 MeV is used throughout the simu- 
lation and the beam kinetic energy of uranium nuclei is 
set to 520 MeV/nucleon if not specifically indicated. In 
the next section, we discuss about the stopping power 
ratio and selection of impact parameter b. In Sec. 3, 
the evolution of baryon and energy densities as well as 
thermalization of central collision systems are studied. In 
Sec. 4, some experimental observables, such as nucleon 
sideward flow and elliptic flow are also investigated. We 
summarize our results in Sec. 5. 



II. STOPPING POWER OF TIP-TIP AND 
BODY-BODY COLLISIONS 

Large stopping power can lead to remarkable pressure 
gradient in the compressed dense matter. It is generally 
also considered to be responsible for transverse collec- 
tive motion [25j], the maximum attainable baryon and 
energy densities as well as thermalization of collision sys- 
tems. Thus, the study of the stopping power in U+U 
collisions may provide important information for under- 
standing the nuclear EoS and collision dynamics. 



A. Selection of impact parameter 

The nuclear stopping power and geometric effects in 
U+U collisions rely strongly on the impact parameter 
b. Considering the conceptual design of the CSR-ETF 
detector [lfjj, two methods are invoked here to estimate 
the impact parameter. The first one is the multiplicity 
of forward neutrons with polar angle 6 < 20° in the lab 
frame which can be covered by a forward neutron wall. 
The other method is to make use of the parameter E rat 
[26| , which is the ratio of the total transverse kinetic en- 
ergy to the total longitudinal one. The particles are also 
required to be within 8 < 20° in the lab frame, while 
the two qualities are calculated within the center of mass 
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FIG. 2: Upper: Forward neutron multiplicity and Lower: 
E rat , as a function of normalized impact parameter b/b max in 
both tip-tip and body-body collisions. 



system (c.m.s.) 



E r 



(1) 



The normalized impact parameter b/b max is used to rep- 
resent centralities of tip-tip and body-body collisions and 
the bmax of the two cases arc quite different from each 
other. As shown in Fig. [2 with either method, obvious 
linear dependence of the normalized impact parameter 
are demonstrated in both tip-tip and body-body near 
central collisions. Then, the two methods can be com- 
bined to determine the impact parameter to identify the 
most central collision events in both tip-tip and body- 
body collisions. 



B. Stopping power ratio definition and evolution 

It is difficult to obtain a universally accepted estimate 
of the nuclear stopping power in heavy ion collisions due 
to a proliferation of definitions of the concept [27]. The 
stopping power ratio R [281 ] is employed to measure the 
degree of stopping and defined as: 



(2) 



, the total nucleon transverse momentum \Pt -\ divided by 
the total absolute value of nucleon longitudinal momen- 
tum \PgA in the c.m.s.. The ratio is wildly used to de- 
scribe the degree of thermalization and nuclear stopping 
by low and intermediate energies heavy ion collisions. 
It's a multi-particle observable on an event-by-event ba- 
sis, which for an isotropic distribution is unity. 

Fig. [3] shows the time and normalized impact param- 
eter dependence of the stopping ratio R for three con- 
ditions: tip-tip, body-body and sphere-sphere collisions. 
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FIG. 3: (Color online) (a)The time evolution of the stopping 
ratio R in tip-tip, body-body and sphere-sphere central colli- 
sions, and (b) the stopping ratio R as a function of b/b max in 
minimum biased collisions. 

When the ratio R reaches the value of 1, full stopping of 
the collision system is considered to be achieved, and the 
momenta is also isotropy, which are not sufficient but nec- 
essary for thermal equilibrium of collision systems [28| . 
For R > 1, it can be explained by preponderance of mo- 
mentum flow perpendicular to the beam direction [2^| . It 
is shown that all of the three conditions can achieve full 
stopping when the stopping ratio R=l, the correspond- 
ing time for body-body and tip-tip central collisions are 
about 15 fm/c and 25 fm/c, respectively. Larger stopping 
ratio and faster evolution to full stopping are observed 
for body-body central collisions than tip-tip and sphere- 
sphere ones at the early stage, which may indicate a more 
violent colliding process for body-body central collisions 
due to the sizable initial transverse overlap region. Al- 
though the stopping ratio of tip-tip central collisions is 
lowest than the other two cases at the early time, it raises 
sharply later and even exceeding one. So, it means that 
longer reaction and passage time can be obtained in tip- 
tip central collisions than body-body and sphere-sphere 
ones, which may indicate the nucleons in tip-tip colli- 
sions can undergo more binary collisions to reach higher 
transverse momentum. 

In Fig. [3Jb), the R of the three conditions are gradu- 
ally decrease with the increase of the normalized impact 
parameter. When b/b max < 0.5, the ratio is always larger 
for tip-tip collisions than the other two cases, while for 
b/bmax > 0.5 all of the three conditions almost have the 
same stopping power ratio. 



FIG. 4: The evolution of (a) baryon and (b) energy densities 
in tip-tip, body-body and Au+Au central collisions. 

to study further about the baryon and energy densities 
evolution in both cases. As the full stopping and de- 
formation effects in U+U collisions, it is believed higher 
local baryon and energy densities system with long du- 
ration can be created, which is considered to be a signif- 
icant condition to study the nuclear EoS at high bayonic 
density region. 

A. The evolution of baryon and energy densities 

The evolution of baryon and energy densities in the 
central zone of tip-tip and body-body as well as Au+Au 
central collisions are illustrated in Fig. [4] 

In Fig. 31 it is observed the maximum attainable 
baryon and energy densities for both tip-tip and body- 
body central collisions are about 3.2 po and 0.8 GeV/fm 3 , 
respectively, while the Au+Au one are about 2.6 po and 
0.6 GeV/fm 3 . Both the baryon and energy densities in 
U+U collisions are higher than the Au+Au one. Once a 
baryon density threshold of p > 2.5 po is required, the 
corresponding duration in tip-tip central collisions ~ 20 
fm/c (from ~ 8 fm/c to ~ 28 fm/c) is longer than ~ 
10 fm/c ( from ~ 8 fm/c to ~ 18 fm/c ) of body-body 
one, which is as predicted. But the peak densities have 
no significant discrepancy between the two cases unlike 
those at the energy region of the Alternating Gradient 
Synchrotron (AGS) [12 ]. which may be attribute to the 
full stopping at the CSR energy. 

B. Thermalization of the U+U collision systems 



III. BARYON, ENERGY DENSITY AND 
THERMAL EQUILIBRIUM 

Considering the discrepancy of stopping power be- 
tween tip-tip and body-body collisions, it is interesting 



As mentioned before, the stopping ratio R = 1 is a 
necessary but not sufficient condition for thermal equi- 
librium of the collision system. In order to approach a 
thermal equilibrium, a long duration of reaction is needed 
for nucleons to undergo sufficient binary collisions. As 
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FIG. 5: The evolution of (a) volume with high density (p > 
2.5po) in tip-tip, body-body and Au+Au central collisions, 



and (b) the scaled mean kinetic energy | < Ek >, within a sions. 
sphere of radius 2/m around the system mass center. 



FIG. 6: Evolution of the multiplicity of the free pion, iV* + 
A, N* + A + tt in (a)tip-tip ,and (b)body-body central colli- 



shown in Fig. Ufa), obvious long duration has been ob- 
tained in both tip-tip and body-body central collisions. 
It is therefore possible thermal equilibrium at the time 
of freeze-out can be achieved. 

The Fig. EJa) is the evolution of the volume with the 
high baryon density(p > 2.5 po) for tip-tip, body-body 
and Au+Au central collisions, respectively. Both tip- 
tip and body-body central collisions have larger volumes 
than Au+Au one at the same beam kinetic energy 520 
MeV/nucleon. Although the maximum volume attain- 
able for body-body central collisions( ~ 220 fm 3 ) is about 
two times larger than tip-tip one(~ 120 fm 3 ), the peak 
volume of tip-tip central collisions lasts a much longer 
time of ~ 10 fm/c (from ~ 15 fm/c to ~ 25 fm/c) and 
much more stable than body-body one. To estimate the 
temperature at the freeze-out time, the scaled mean ki- 
netic energy of all hadrons in a sphere of radius 2/m 
around the system mass center is calculated as | < Ek > 
[22| . which is utilized to reflect the thermalization tem- 
perature T of the collision system approximately. As 
illustrated in Fig. [5] (b), both tip-tip and body-body 
central collisions show a flat region about 75 MeV and 
the corresponding time range are about 10 fm/c to 28 
fm/c and 10 fm/c to 18 fm/c, respectively. Considering 
the time range of the flat region in Fig. [5] (b) associat- 
ing with the corresponding range in Fig. [5] (a) and also 
looking back to Fig. QJ we obtain a large volume of hot, 
dense nuclear matter in both tip-tip and body-body cen- 
tral collisions. Consequently, the extreme circumstance 
of sufficiently high temperature and densit y fo r a signif- 
icant large volume and long duration (l2l l22l| has been 
formed in tip-tip and body-body central collisions, which 
can provide a good opportunity to study the nuclear EoS 
as well as particles in medium properties, especially for 
tip-tip case. 

The time of freeze out should be cautiously determined 



for estimating the thermalization temperature of colli- 
sion system. In Fig. [6] the multiplicity evolution of free 
pion which are not bounded in baryon resonances and 
pion still bounded inside the excited baryon resonances 
(A, N*) (unborn pion) are displayed. At the Lanzhou 
CSR energy region (520 MeV/nucleon), the production 
and destruction of the A resonances are mainly through 
AW ^ NA and A — > Ntt reactions in which the A decay 
rate is always higher than that of the formation of this 
resonance and the production of pion is predominated by 
the decay of the A resonances (A — > Ntt) [3(J. The total 
multiplicity of pion, A and N* approaches a saturated 
level after a period of evolution, indicating the freeze- 
out time about t=28 fm/c and t=18 fm/c for tip-tip 
and body-body central collisions, respectively. The larger 
maximum attainable total multiplicity of pion, A and N* 
and freeze out earlier indicates a existent of faster evo- 
lution and more violently reaction process for the body- 
body central collisions than tip-tip case consisting with 
the discussing before. 

The corresponding temperature about 75 MeV at 
freeze-out time can be extracted from the Fig. [S](b), for 
both tip-tip and body-body central collisions. To further 
confirm this estimation, both the energy spectrum of the 
nucleon and negative-charged pion are studied within the 
polar angle range of 90° ± 10° in the c.m.s.. The thermo- 
dynamic model [3lj predicts that the energy spectra will 
be represented by a temperature T which characterizes a 
Maxwell-Boltzmann gas 
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, where P and E are the particle momentum and total 
energy in the c.m.s.. Both the energy spectra and the 
Boltzmann fit results are shown in Fig. [7] The inverse 
slope (e.g. temperature T) of the nucleons in tip-tip and 
body-body central collisions are about 73 MeV and 70 
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FIG. 7: (a) Nucleon, and (b) negative-charged pion energy 
spectrum at 90° ± 10° in the c.m.s. together with a Maxwell- 
Boltzmann fit for both tip-tip and body-body central colli- 
sions. The nucleon fit temperature for tip-tip and body-body 
are about 73 MeV and 70 MeV, respectively and that of pion 
are about 56 MeV and 52 MeV, respectively. 



MeV, respectively, which are in good agreement with the 
temperature extracted from the Fig. [5jb) at the freeze 
out time. The spectra of negative-charged pion show a 
different lower temperature than that of nucleon which 
may be explained by considering an equilibrated N and 
A system at thermal freeze out and taking into account 
the kinematics of A decay (32j . The nucleon temperature 
closely reflects the freeze-out temperature of tip-tip and 
body-body central collisions. 

In conclusion, thcrmalization (or near thermalization) 
of the collision system corresponding a freeze-out tem- 
perature about 75 MeV is likely to be achieved in both 
tip-tip and body-body central collisions. However, it's 
also possible that the collision system is still in a non- 
equilibrium transport process on its path towards kinetic 
equilibration (30j . 



IV. THE COLLECTIVE FLOW OF U+U 
COLLISIONS 

Stopping of nuclei in heavy ion collision can lead 
to pressure gradient along different directions, result 



ing in collective motion as spectators bounce-off [34 1 
and participants squeeze-out effects (35|. Since last 
two decades, at Bevalac/LBNL and SIS/GSI energies 
the so-called "collective flow" analysis has been estab- 
lished [H, [H, HE H H3 to study the collective mo- 
tion of the products in heavy ion collisions. The collec- 
tive flow resulting from bounce-off and squeeze-out ef- 
fects, which can be explained well by the hydrodynamics 
model [34 , 1 3 81 ] , and also be in good agreement with the 
experimental data has been observed [3^, H(| ■ Because of 
the large deformation of the uranium nuclei, a novel col- 



FIG. 8: The mean transverse momentum per nucleon pro- 
jected into the reaction plane, < p x /A >, as a function 
of c.m.s. normalized rapidity is illustrated for tip-tip and 
body-body collisions. With normalized impact parameter 
cutoff: (a)b/b m ax <= 0.5 (h)b/b max > 0.5. 



lective motion is expected [l2|, to be used to extract the 
medium properties and nuclear matter EoS information. 

To perform flow analysis, it is necessary to construct a 
imaginary reaction plane defined by direction of the beam 
(z) and the impact parameter vector b (43l. |45|. I4r3| . In our 
simulation, the x — z plane is just defined as the reaction 
plane with the beam direction along z positive direction 
and the impact parameter vector b along x positive direc- 
tion. In last two decades, there are mainly two methods 
to study the collective flow at the low and intermediate 
energies. One is the sphericity method [H, [H, HJ, |42j 
which yields the flow angle relative to the beam axis of 
the major axis of the best-fit kinetic energy ellipsoid, and 
the other is to employ the mean transverse momentum 
per nucleon projected into the reaction plane, < p x /A >, 
to perform nucleon sideward flow analysis [43l . |44| which 
reflects the spectator bounce-off effects in the reaction 
plane. In recent years, it is usual to use an anisotropic 
transverse flow analysis method. With a Fourier expan- 
sion [47], EH of the particle azimuthal angle <p distribu- 
tion with respect to the reaction plane, different har- 
monic coefficients can be extracted, among which the 
first harmonic coefficient v\, called directed flow (simi- 
lar to sideward flow) and the second harmonic coefficient 
t>2, called elliptic flow are mostly interested. The ellip- 
tic flow reflects the anisotropy of emission particles in 
the plane perpendicular to the reaction plane while the 
directed flow describes the anisotropy in reaction plane. 
The Fourier expansion can be expressed as 
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2v n cos{n(f>) 



(4) 



Fig. [5] shows nucleon sideward flow, < p x /A >, for 
both tip-tip and body-body minimum biased collisions 
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used quality to describe the nucleon sideward flow quan- 
titatively defined as 



b/b m 



FIG. 9: (a)The nucleon flow parameter F and (b)the c.m.s. 
mid-rapidity ( —0.5 < y° < 0.5 ), nucleon elliptic flow V2 of 
three collision conditions as a function of normalized impact 
parameter b/b max with soft EoS. 



as a function of normalized rapidity, y^ = Y crn /y cm , in 
which Y cm represents the particle rapidity in c.m.s. and 
y cm is the rapidity of the system mass center. To ex- 
tract the nuclear EoS information and also demonstrate 
the discrepancies of the nucleon sideward flow in tip-tip 
and body-body collisions, the cascade events [49(, which 
neglect the mean field and pauli blocking effects are em- 
ployed here to compare with the soft EoS case. In Fig. 
[HIa),(b), with a soft EoS, it is noted that either tip-tip 
or body-body collisions show a spectator bounce-off ef- 
fect revealing an obvious " S" shape [H, E3| at the mid- 
rapidity region of— 0.5 < y° < 0.5, while the cascade one 
appear a almost vanishing nucleon sideward flow. It can 
be understand by the nucleon sideward flow is related 
to the mean field, which is mainly responsible for the 
pressure gradient of the stopping nuclei, while the mean 
field has a strong dependence of the nuclear EoS. There- 
fore, the nucleon sideward flow is thought to be a good 
indirect probe to extract the nuclear EoS information, 
especially tip-tip case for its largely remarkable sideward 
flow. A cutoff on normalized impact parameter is also 
applied to explore the impact parameter dependence of 
nucleon sideward flow. As shown in Fig. [Hb), when 
b/b max > 0.5 the curves of soft EoS and cascade are al- 
most superposed with each other, while for b/b ma x < 0.5 
large discrepancy is observed. The situation is quite sim- 
ilar to Fig. [3] (b), almost the same stopping power for 
b/b m ax > 0.5 and large discrepancy for b/b max < 0.5 in 
tip-tip and body-body minimum biased collisions, which 
means there exists a correlation between nuclear stopping 
power and sideward flow (33j . 

The normalized impact parameter dependence of the 
collective flow of nucleon is further studied, by analyzing 
the "flow parameter" F [4^| and also elliptic flow t>2 for 
both tip-tip and body-body as well as Au+Au minimum 
biased collisions. The flow parameter F is a customarily 



F 



d < p x /A > 



dy(°) 



y(0)=0 



(5) 



the slope of the mean transverse momentum per nucleon 
projected into the reaction plane at z/°) = 0. 

In Fig. [9]Ja), with b/b max > 0.5, the nucleon flow pa- 
rameter F of tip-tip and body-body collisions are with 
similar value. This similarity, along with the almost same 
stopping ratio R in Fig. [3]Jb), indicates a existence of sim- 
ilar pressure gradient effects on nucleon sideward flow in 
the two collision orientations. While for b/b max < 0.5, 
the flow parameter F of tip-tip collisions is nearly 3 times 
larger than that of body-body case. Even the sideward 
flow of Au+Au collisions is larger than the body-body 
one. It is further confirmed the tip-tip nucleon sideward 
flow is a more sensitive probe to extract the information 
of nuclear EoS than that of body-body one. The promi- 
nence high of the nucleon sideward flow in tip-tip colli- 
sions may be resulted from the stronger pressure gradient 
between the participants and spectators in the reaction 
plane than body-body one, due to the largely deformed 
nuclei. 

The normalized impact parameter dependent of nu- 
cleon elliptic flow V2 at the mid-rapidity region ( —0.5 < 
y° < 0.5 ) is displayed in Fig. \9lb). A significant neg- 
ative elliptic flow V2 at this energy region is consistent 
with the excitation function of the elliptic flow studied 
before [501 ] . An largest negative V2 about —12% in body- 
body central collisions is observed which reflects the large 
geometric and squeeze-out effects in the collisions. While 
for tip-tip and Au+Au ones the maximum negative V2 are 
obtained at mid-centrality. Since both high baryon, en- 
ergy densities and large elliptic flow effects, which reflects 
an early EoS of the hot dense compression nuclear matter 
[r^ . are available in body-body central collisions. Thus 
the body-body nucleon elliptic flow can also be taken as 
a sensitive probe of nuclear EoS. The novel behaviors 
of nucleon collective flow in tip-tip and body-body col- 
lisions are mainly attributed to the large deformation of 
the uranium nuclei. 



V. SUMMARY 

In summary, the CSR-ETF at Lanzhou provide a good 
opportunity to systematically study the nuclear EoS at 
the high net-baryon density region of nuclear matter 
phase diagram. Due to the novel stopping effects in 
largely deformed U+U collisions, the simulation based on 
ART1.0 demonstrates that full stopping can be achieved 
and also a bulk of hot, high densities nuclear matter with 
large volume and long duration have been formed in both 
tip-tip and body-body collisions. Large nucleon sideward 
flow in tip-tip collisions and the significant negative nu- 
cleon elliptic flow in body-body central collisions can pro- 
vide a sensitive probe to extract nuclear EoS information. 
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Thus the extreme circumstance and the novel collective 
flow in both tip-tip and body-body collisions can provide 
a good condition and sensitive probe to study the nu- 
clear EoS, respectively. More experimental observables 
of U+U collision dynamics should be further studied, due 
to the geometric effects. 
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